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Abstract 



We present a QCD analysis of the strange and charm contributions to the neutrino 
deep inelastic structure function xF 3 . We show that next-to-leading order effects, which 
are relatively important for F2, play a lesser role in the case of xF%. The neutrino- 
antineutrino difference xF% — xF^ provides a new determination of the strange density, 
which exhibits some advantages with respect to other traditional methods. 
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Two methods are traditionally used to extract the strange-sea density from deep 
inelastic scattering (DIS) data: the first consists in studying charm production in 
charged- current neutrino DIS (the characteristic signature of this process being the 
presence of dimuons in the final state); the second is to subtract the Fi structure 
functions measured in neutrino and muon DIS, thus selecting the strange contribution. 

Until last year, these two determinations, based on the NMC /iDIS data JIJ and 
on the CCFR z/DIS data 0, [| available at that time, seemed to yield contradictory 
results for s(x). This discrepancy, which has strongly challenged the attempts at global 
parton fits [£|, H, had been actually predicted || 0], and is simply explained || |j| by 
the observation that the dimuon and v — [i determinations of s[x) actually measure 
different quantities, related to - but not coincident with - the strange density. This 
is due to the relevance of quark-mass effects Q and longitudinal contributions, in the 
region of small and moderate Q 2 values (of order of 10-30 GeV 2 ) fj 

A recent QCD analysis [H| of the new CCFR dimuon data [[l l[ has shown that, 
when all important physical effects are taken into account, the different measurements 
converge - as they should - towards a unique result for s(x). Incidentally, the strange- 
quark distribution emerging from all data is well reproduced by a "traditional" fit, 
such as, for instance, MRS(A) |T2|], and does not support a nearly SU(3) symmetric 
fit, such as CTEQ1 §. 

Both the dimuon and the v-[i extractions of s(x) present problems and subtleties 
(for a detailed discussion see |9|, [HI). In particular, the dimuon determination implies, 
for experimental reasons, an acceptance-dependent separation of the t- and ^-channel 
diagrams that constitute the 1^-gluon fusion QCD process. On the other hand, the 

v — fi result is affected by large uncertainties due to the very unsafe procedure of 

2 Remember that in charged-current neutrino DIS strange and charm excitations are inseparable. 
3 For instance, the average Q 2 value of the CCFR data is {Q 2 ) ~ 22 GeV 2 . 
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subtracting data from two different experiments. 

However, the idea of obtaining s(x) by an appropriate combination of structure 
functions can be further exploited. There is in fact another way to isolate s(x) from 
DIS structure functions, which makes use of the third z/DIS structure function, F3. In 



the parton model these are [13 



xF£ N (x) = xV(x) - 2xc(x) + 2xs{x) , (1) 
xF* N (x) = xV(x) + 2xc{x) - 2xs{x) , (2) 

where V(x) is the valence distribution and N denotes an isoscalar nucleon (hereafter 
we shall drop the suffix N from our formulas). Therefore, the v—v difference effectively 
measures the strange density, since the charm contribution is very small, at least in 
the kinematical region investigated by the present experiments (we assume s = s and 
c = c): 

xF% N {x) - iff (1) = Ax [s{x) - c{x)\ . (3) 

Needless to say, moving from the parton model to leading-order QCD, the quark dis- 
tributions acquire a Q 2 dependence governed by the Altarelli-Parisi equations. 

The use of eq. (§) to extract xs(x) has the immediate advantage, over the v — fi 
method, of combining data from the same experiment, thus with no relative- normalization 
problems. However, in practice there are at least two shortcomings. First of all, xF 3 is 
a small, valence-dominated, quantity, more difficult to measure than F2. Secondly, the 
statistics for xF^ — xF^ is much lower than that for F 2 : the measurement of xF£ — xF% 
requires neutrino and antineutrino data separately, whereas they can be combined for 
F2, which is the same in uN and vN deep inelastic scattering. In general, the number 
of //-induced events is considerably smaller than that of ^-induced events: for instance, 
the ratio is about 1 to 5 in the CCFR experiment 
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The CCFR/NuTeV Collaboration is at present working on the extraction of xs(x) 
from xF^ N — xF^ , and data will be available in the near future [|14|j . 

Previous studies on the determination of the strange density at moderate Q 2 from 
neutrino DIS have taught us the importance of quark-mass corrections and current non- 
conservation effects, which manifest themselves through the order-a s vector-boson- 
gluon fusion diagrams. These represent the dominant contribution near the heavy- 
quark threshold, in particular at small x. It is then natural to go beyond the leading 
order also in the analysis of the extraction of xs(x) from xF%. The Next -to-Leading 
Order corrections are known to affect mostly the longitudinal component of structure 
functions 0: thus we expect the NLO effects to be smaller in xF 3 than in F 2 , because 
xF% is a purely transverse structure function. However, only an explicit computation 
can give us precise information on the charm-strange content of xF 3 . 

In the following we shall present a QCD calculation of xF^—xF^ at order a s , taking 
into account the contribution of the gluon-fusion diagrams. The possibility of a safe, 
unambiguous, extraction of the strange density from xF^ N — xF^ N will be explored 
and discussed. 

At order a s the main contribution to the cs component of F% is given by the W— 
gluon fusion (GF) term 4 , which, for the strange-charm sector, reads 

Here a = 1 + (m 2 + m' 2 ) / Q 2 and, for neutrino scattering, m is the mass of the charmed 
quark, m! is the mass of the strange antiquark. The Wilson coefficient C3 represents 
the W + g —>■ cs cross section difference 0^ — 0^ (L and R standing for left- and right- 



transverse W, respectively); it has the explicit form [15], [16] 

( 2/2 
Ttl — Ttl 

C 3 (z,Q 2 ) = 2 p + /3_ ^ 2z(l - z) 



Q 2 



^The 0(a s ) VF-quark fusion diagrams are a negligible correction. 



2 Q 2 ^ Q 4 " J ~ v ""' ' 

+ \m <r-¥ m'] L(m',m)} , (5) 



where 



and 



1 + + /?+/?- 

L(m,m) = log 2 , 2 ^ . (7) 

v ' ' 1 i m z —m' z zoo v ' 

1 H q5— — - P+P- 
If the two quark masses are non-zero, the Wilson coefficient C3 is free from singu- 
larities. In the limit ml — *> 0, i.e. treating the strange quark as massless, L(m,m') and 
L(m',m) behave as 

g 

L(m,m') — ► log — - (8) 
m 2 

(g — Tfl 2 ) 2 

L(m',m) -> log v ; ,/ (9) 
s m' z 

where s = Q 2 (l — z)/z. In the same limit, the factor multiplying L(m',m) becomes 
V g (£), namely the usual g — > qq splitting function V g expressed in the rescaled variable 
£ = x(l + m 2 /Q 2 ). While L(m,m') is regular, L(m',m) has a collinear singularity. 
This is subtracted out by setting JL7]] 

L(m',m) = log (g T^ 2)2 , (10) 

where the scale /z 2 is customarily taken to be equal to the factorization scale (i.e. to 
the scale that separates the perturbative part from the non-perturbative one in the 
DIS QCD diagrams). 

Let us now come to the total cs contribution to F 3 . We start from the QCD 
factorization formula, which formally reads (<g means convolution) 

F(x, Q 2 ) = £ Mz, /i 2 ) <g> C(~ /i 2 , Q 2 ) , (11) 



where the sum is made over all parton species. 

At leading order only quark and antiquark contribute and the corresponding Wilson 
coefficients are delta functions of the slow-rescaling variable £ = x(l+m 2 /Q 2 ). At next- 
to-leading order the contribution is given by eq. (Q) with the above-defined subtraction. 
Thus we get 

'a,\ f 1 dz , n. ( 3c 



s(Z^ 2 )-c(x^ 2 )} + (^-) [ -g(z^ 2 )C 3 (-,Q 2 ) . (12) 

J V 71 J J ax Z \Z J 



Here C3 stands for the subtracted Wilson coefficient, i.e. for C3 with the replacement 
(|T0|). Notice that the quark excitation (QE) term, s — c, is taken at the factorization 
scale /i 2 . The Q 2 evolution (at least the dominant part of it, due to g — > qq splitting) 
is already contained in the gluon-fusion Wilson coefficient. Both taking the quark 
densities in eq. (0) at the physical scale Q 2 and using C 3 instead of C3 in eq. 
would represent a double counting. A most often used approximation consists in setting 

Fr(x,Q 2 ) = 2[s(Z,Q 2 )-c(x,Q 2 )] , (13) 

which means combining massless QCD, for the Q 2 evolution, with slow rescaling, to 
account for quark mass effects. In the following we shall check the goodness of this 
slow rescaling procedure for F3. 

At large momentum transfer (Q 2 ^> m 2 , say Q 2 > 30 GeV 2 ) the charmed quark 
should also be treated as a massless parton. This means that L(m, m') too becomes 
large and should undergo a subtraction similar to that performed by eq. (0). Asymp- 
totically, the massless QCD formulas are of course regained: in the limit m, m' — > 
the contribution of the gluon-fusion diagram to F3 vanishes, i.e. F^q F 0. 

So much for the neutrino DIS. In the case of antineutrino scattering the formulas 
written above must be modified by exchanging m and ml: m becomes the strange mass, 



m' the charm mass. The GF contribution changes sign and so does the QE term. The 
difference of the whole structure functions, F£ — F£, is simply twice the cs component 
for neutrino, 2F^ ,CS , since the valence part cancels out (notice that the cs component 
disappears in the sum Fg + Fg and hence is irrelevant for the Gross-Llewellyn Smith 
sum rule). 

Let us now present the results of our calculations. The v — v difference of xF% 
structure functions, xF% —xFg, has been evaluated by using eq. flT2| ) and the MRS(A) 



fit [12] for the strange, charm and gluon densities at the factorization scale fi 2 . The 
choice of /i 2 is a delicate issue. As we shall see, one of the advantages of working 
with xF£ s is that its dependence on fi 2 turns out to be quite small (smaller than the 
/^-dependence of F™). In any case, /i 2 = m 2 , where m c is the charm mass, seems to 
be a natural choice near (or not much above) threshold. At large Q 2 it is reasonable 
to take a factorization scale of order Q 2 , rather than m 2 . The prescription of Ref. [17], 
for instance, gives /i 2 ~ m 2 at small Q 2 and /i 2 ~ Q 2 /2 at Q 2 3> m 2 . 

In Fig. 1 we present our results for xF£ — xF£ just above threshold (Q 2 = 10 GeV 2 
and 25 GeV 2 ), with /i 2 = m 2 (we use m 2 = 2.7 GeV 2 ). In this case c(x,fi 2 ) obviously 
vanishes. The total result, eq. (0), the quark excitation term, the unsubtracted gluon 
fusion contribution, eq. (|), and the slow-rescaling prediction, eq. (|T5|), are plotted in 
the figure. An interesting feature is clearly visible: the complete NLO result for xF% — 
xF£ nearly coincides with the slow-rescaling expectation, whereas it differs sensibly 
from the unsubtracted GF result, especially for x > 0.01. Thus, the slow rescaling 
mechanism, whose application to the longitudinal+transverse structure function F2 at 
small Q 2 is rather unsafe (see 0, 0), turns out to be an excellent approximation when 
dealing with the charm-strange contribution to F3. This is due to the fact that the 
main next-to-leading order effects that slow rescaling mimicks too crudely are related 
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to the longitudinal component of structure functions, which is absent in F3. Notice 
also that the QE component, which is the only term containing the strange density 
that experiments aim to extract, is comparable in magnitude to - and actually not 
much different from - the complete result. 

We checked the dependence of the results on the factorization scale fi 2 . This de- 
pendence is shown in Fig. 2 for two values of x (0.01 and 0.1) and for Q 2 = 25 GeV 2 
(close to the average Q 2 value of the CCFR experiment). It is reassuring to see that, 
not only the total result for xF% — xF%, but each single term (and in particular the 
QE term) has a very mild dependence on /i 2 . Moreover, the (small) fj 2 - dependence of 
the QE term is approximately the same as that of the full v — v difference (the solid 
and dotted curves in the left part of Fig. 2 are parallel). For comparison, in the right 
part of Fig. 2 we present also the scale dependence of F%' cs , which appears to be more 
dramatic, in particular in the QE component (for more details we refer the reader to a 



forthcoming paper fll8[). Since it is the latter component that contains the strange and 
charm densities, it is evident that the extraction of xs(x) from xF£ — xF£ is affected by 
a factorization scale uncertainty much smaller than that occurring in a measurement 
based on F^' 08 , such as the dimuon measurement. 

At large Q 2 both strange and charm can be considered as massless partons and one 
expects to regain the results of massless QCD: in particular, the subtracted GF term 
should vanish. This is clearly visible in Fig. 3, where one can see that the full result 
coincides asymptotically both with the QE term and with the slow-rescaling prediction 
(this means that F^'q f - represented by the dashed curve - is exactly cancelled by the 
subtraction term, so that F^'q F = 0). 

Leaving aside the lack of statistics that may make the determination of the strange 
density from xF% — xF% difficult in practice, it is clear that this method has some 
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indisputable advantages: i) it is not plagued by relative-normalization errors; it) it 
is not affected by the ambiguities inherent in other methods, such as the spurious 
separation of t- and ^-channel diagrams that occurs in the dimuon separation (see 
H [UJ); Hi) large longitudinal contributions arising from the non-conservation of weak 
currents are obviously absent; iv) charm mass effects are very well accounted for by 
the slow-rescaling prescription, making the analysis of data and the extraction of the 
strange density particularly simple (all the next-to-leading order QCD machinery can 
be safely avoided); v) the dependence on the factorization scale of the full 0(a s ) result 
and of the quark excitation term, which contains the strange density to be extracted, 
is rather mild and does not represent a worrisome source of uncertainty 

The main conclusion of our study is that a precision measurement of xF% — xF% 
may well represent a new, interesting source of information on the strange quark distri- 
bution: being cleaner and more direct than other determinations, it is certainly worth 
exploiting. 

Useful discussions with A. Bazarko and E. Predazzi are gratefully acknowledged. 
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Figure Captions 



Fig. 1 The difference xFg - xF% at Q 2 = 10 GeV 2 and 25 GeV 2 . The solid curve is the 
complete result, eq. dl2|) . The dotted and dashed curves are the quark excitation 
(QE) and the unsubtracted gluon fusion (GF) contributions, respectively. The 
dot-dashed curve is the slow-rescaling expectation, eq. fljjl) . The factorization 
scale is /i 2 = m 2 . 

Fig. 2 The dependence of xF% — xF% and of F^' 08 on the factorization scale /i 2 . The 
meaning of the curves is the same as in Fig. 1. 

Fig. 3 Same as Fig. 1 at Q 2 = 100 GeV 2 and 1000 GeV 2 , with /i 2 = Q 2 /2. Note that in 
the bottom window the solid, dotted and dot-dashed curves coincide. 
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This figure "figl-2.png" is available in "png" format from: 
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